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I. INTRODUCTION 

This paper concerns antiferromagnets with isotropic ex- 
change couplings, i.e. the Hamiltonian has form 

Ti-cx — —-^ JjjSj ■ Sj. (1.1) 

(Such sums will always run to N, the number of spins.) When 
the interactions in il.li are frustrated, the classical ground 
states often exhibit continuous nontrivial degeneracies — be- 
yond those induced by the global spin symmetry. In a highly 
frustrated antiferromagnet', the number of such degrees of 
freedom is extensive; in other cases, such as face-centered cu- 
bic (fee) antiferromagnets, it is a finite-dimensional manifold. 
In either case, such degeneracies are lifted by thermal fluc- 
tuations, quantum fluctuations, or quenched fluctuations (due 
to dilution). The selection of a ground state due to fluctua- 
tions was called^ (by a slight abuse of terminology) ordering 
by disorder.^ The purpose of this paper is to provide effec- 
tive Hamiltonians, in convenient form, which represent the 
degeneracy-breaking selection (free) energy. 

A ground state selection calculation for a continuous spin 
system was performed first by Shender, who showed that 
in the two-sublattice garnet Mn3Cr2Ge30i2 quantum fluc- 
tuations favor a coUinear spin configurationii The selection 
causes an effective "anisotropy" of one sublattice relative to 
the other, which opens a gap in some soft spin wave modes ob- 
servable in neutron scattering.''-^ Others have calculated a sim- 
ilar effect for the type I FCC antifeiTomagnet 7-MnpFei_p^ '^ '' 
. Rastelli and co-workers investigated thermal and quantum 
fluctuations in the frustrated rhombohedral lattice (describ- 
ing f3 — O2), finding again that quantum fluctuations select 
collinear ground states'". The collinearity bias was deci- 
sive between different possible fee orderings in Cu nuclear 
spins''"'^. Quantum fluctuations were also found to decide 
the ground state in Sr2Cu304Cl2^. 

Whereas thermal or quantum fluctuations select collinear 
states"' '-' quenched fluctuations due to dilution select non- 
collinear states^J^, in cases where the ground state manifold 
is finite dimensional. (In highly frustrated magnets, meaning 



the degrees of freedom are extensive, the many zero-energy 
excitations "screen" defects and dilution does not necessar- 
ily favor a global noncollinear stat o'^i'^ .) The competing 
selection effects can give a rich phase diagram as a func- 
tion of temperature and dilution, as simulated in the J1-J2 
square lattice antiferromagne t ' 3, 1 8,20 _ jjj particular, mak- 
ing use of the effective Hamiltonians presented in this pa- 
per, we studied the fee type III antifeiTomagnet with dilu- 
tion and Dzyaloshinskii-Moriya anisotropic exchangsiS; this 
represents the antiferromagnetic dilute magnetic semiconduc- 
tor Cdi-pMnpTe for x — > 1, a material realized using 
molecular-beam epitax y^ ''^^ . 

For all three kinds of selection - thermal, quantum, and 
dilution - the selection is conveniently expressed by adding 
biquadratic exchange terms to the Hamiltonian, 

ij 

Such a form is quite familiar for selection by quantum fluc- 
tuations in the "independent sublattices" case, i.e. the an- 
tifeiTomagnet consists of two (pr more) sublattices with a 
strong, unfrustrated intra- lattice intereaction Ji, and a weak, 
frustrating inter-lattice interaction j^iSiiliSi: Eq. ( 11.21 ) fol- 
lows directly from perturbation in Ji / J2, with inter-sublattice 
Kij — 0{J1/ J2). However, there are many other cases in 
which the degeneracy freedom requires correlated rotations in 
all the sublattices. This paper lays out the framework in which 
( 11.2b can be justified (approximately) in that broader class of 
frustrated systems. 

In Sec. [n] we present a microscopic calculation of the free 
energy in the presence of quenched disorder and thermal fluc- 
tuations, performed to lowest nontrivial order in spin devia- 
tions. The combined selection effects are extracted in a new 
way by performing a constrained integration over the fluctua- 
tions in reciprocal space. Then, in Sec. lIIBI we express the re- 
sults in the phenomenological form of an effective biquadratic 
exchange term. In Sec. |III1 quantum fluctuations are shown 
to give a similar form, using an approximation of the spin- 
wave approximation. Appendices |A]and|B]give, respectively, 
an alternate path to the dilution effective Hamiltonian, and a 
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calculation of magnetoelastic effects (which reduce to a sim- 
ilar quartic effective Hamiltonian, and thus give selection ef- 
fects similar to those of thermal or quantum fluctuations). The 
Discussion (Sec. IIVI ) proposes applications for the effective 
Hamiltonian in working out phase diagrams in simulations of 
magnetic systems. 

II. THERMAL FLUCTUATIONS AND DILUTION 

In this section we show that thermal fluctuations favor 
collinear states while dilution favors the least collinear states. 
This competition is a general feature of all uniformly frus- 
trated vector spin systems studied so far and has a simple 
explanation^ '^. It is familiar that spins in a vector antifer- 
romagnet will orient transverse to an applied magnetic field, 
because they can thereby gain energy by relaxing towards it. 
The effective local fields generated by thermal fluctuations are 
transverse to the fluctuating spins, and the free energy will be 
lower when other spins are transverse to these effective fields. 
The coupling, and the entropy, is therefore maximized when 
all spins are collinear*. 

By contrast, dilution (the removal of spins) generates effec- 
tive fields parallel to the removed spins (or other spins of the 
same sublattice). In this case the energy is minimized for the 
least collinear states, because most spins are again transverse 
to the effective field acting on them. The above qualitative 
arguments will be confirmed by perturbative spin-mode cal- 
culations. 

In the rest of this section, we first set up the artificial form 
of dilution; then, representing a delta function by auxiliary 
variables, we sidestep the awkwardnesses expanding about a 
general nonperiodic, noncollinear classical ground state (Sub- 
sec. UlAl. This result can be put in the form of a biquadratic 
effective Hamiltonian (Sec. lIIBl i. 

A. Formal Derivation of Selection Free Energy 

For this section, the spins in il.lj are taken to be classical 
unit vectors with n components. We compute the free energy 
by assuming the fluctuations about the ground state are small, 

S, -SW, + cr, , S(°\;-cr, = , (2.1) 

and \(t\ ^ 1. Further, we will consider "modulation" disor- 
der where the exchange between sites i and j is increased or 
decreased by the independent random variables e^. The actual 
couplings are taken to be 

J--i = £.)(!- 6,) , (2.2) 

with (e,;) = 0. This approximation^"' lets us control the size of 
the deviations cr, by taking the variance of the random mod- 
ulation to be a small parameter: (e^^) = (e^) <C 1. We will 
take J', — Jf,"'° for modulation disorder, where Jf,"'^'^ is the 
energy of a bond in the pure system. 

True site dilution means J^J^' — Jf^'^ili'riiy, where the ran- 
dom variable rji is 1 (site occupied) with probability p and zero 



with probability zero. (See Appendix lAl) This takes the form 
of Eq. ( 12. 2t , with (e^) — 0, if we take = 1 — rji /p. This 
then forces Jy — p^Jf^^'^. 

For the "independent lattices" case e.g. the type II bcc'*-, 
the type II FCC ordering, the frustrated square lattice modeU^, 
or a body-centered tetragonal antiferromagnet^^, the ratio 
•h/ J2 can be a small parameter for achieving this same con- 
trol in the limit of large | J2I, because the noncollinear part 
of the effective field of a missing spin becomes small com- 
pared to the collinear part. In AppendixlAlit is shown that the 
model of ( 12.21 ) and true dilution are equivalent to lowest order 
in the spin fluctuations. Monte Carlo results for true dilution 
are consistent with the results of this theory for p < lfi^i2£ 

The ground state condition is that sf^ — h|"-' / jh^"'' | for all 
i, i.e. every spin is aligned with its "local field" 

3 

We assume Ih-^"*! = ho, the same at every site (subsequent 
manipulations depend on this.) The assumption is valid on 
most cases that the sites are symmetry-equivalent, even in 
highly frustrated systems (e.g. pyrochlore lattice) where a 
generic ground state is not at all periodic. Substituting ( 12.31 ) 
in dl.lt shows Eq = ho/2, so Hq must be the same in every 
ground state, too. 

Expanding in small deviations, the Hamiltonian dl.ll) be- 
comes a quadratic form 

i i 

+ i^A(z,j>,.<T,+0(e3i)2.4) 

Here E^ is the ground state energy per spin of the pure system, 
and 

h, = ^J,,,ejSf . (2.5) 
3 

Also, 

= —Ji,j + hoSij . (2.6) 

The first term on the right hand side of ( |2.4| i is indepen- 
dent of cr j and vanishes when the configurational averages are 
taken. We ignore the last term since it is higher order in e. 

1. Evaluation via auxiliary representation of constraints 

We now have a Hamiltonian which is purely quadratic and 
apparently trivial. However, as always, each spin really has 
only 71—1 degrees of freedom since its length is fixed. When 
n > 2, there is an arbitrariness in the choice of basis vec- 
tors for the local transverse subspace, a gauge freedom. It is 
cumbersome to fix this gauge by a choice of local frames^-. 
Instead, following Ref. 2 we shall implement the constraints 
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so as to avoid introducing an explicit local frame for each spin, 
using integral identities to maintain a manifestly rotation- 
invariant (gauge-invariant) form. 

For convenience, a matrix notation will be used where the 
matrix indices specify both position (z) and vector component 
in spin space (//). The unit length constraint becomes 



P • (T = 0. 



(2.7) 



Here P is the nonsquare {N x Nn) matrix whose components 
are 



Since S -^"* is a unit vector, 

ppT 



(2.8) 



(2.9) 



where I is the identity matrix, and we (sometimes) add a 
subscript to clarify a matrix's dimension. (However, notice 
P-^P ^ T-NnxNn-) The definition of Aij [Eq. (I2.6I) 1 becomes 



hnl. 



(2.10) 



In ( I2.IOI 1. the matrices A and J are N x N, but from here till 
the end of the section they are extended to be Nn x Nn matri- 
ces (by taking the direct product with Inxn, acting on the spin 
indices). Also, h [from ( 12.51) 1 becomes an A^n-component 
vector in this notation. 

Our goal, the partition function for Hamiltonian ( 12.41 ). is 



Z = J [Y[da.,^][llSiF ■ a) 



cxp ( - i/Jcr'^Acr + ph^a 



(2.11) 

By representing the delta function of the constraints in terms 
of a functional integral of exp(itPcr) over auxiliary variables 
{ti} and completing the square in the resulting exponential, 
the dcTii^i integral can be done^, yielding 



Z = Zi 



exp(- i(h^ + z/3 



-l.T 



(2.12) 



Here 



Zi 



[n 



dcr, 



exp 



1 



-/3cr^ Acr 



(2.13) 



To make further progress, let us define the (N x N) matrix 
of the quadratic coefficients in ( |2.12| i: 



B 



PA^P^, 



(2.14) 



All our results will be expressed in terms of B. Since B is 
nonsingular, we can evaluate ( |2.12| i using a second comple- 
tion of squares^i. 



Thus, 

where 

^2 = 



Z = Z\ Z2 e 



-0AF 



(2.15) 



[n 



dtj 
2^ 



[27r/3]-(detB)-5, 



(2.16) 



and 



AF = -ih^^A-i - A-ip^B-ipA-i)h. (2.17) 



2. Free energy result and disorder average 



Thus the exact free energy is 



1 



F = -p-' A^(l + -) ln(27r) + iV(l - -) ln/3- -ln(det A) 



+ /3-iiln(detB) + Ai^-h^-S('') 



Fo + Ft + F^ii + F] 



REF 



scl ? 



(2.18) 



where we have reintroduced the linear term in h from (12.41 ). In 
this equation, the term in brackets (Fq) is independent of the 
ground state, hence irrelevant for state selection (though pos- 
sibly important in comparing to the free energy of competing 
orders). The remaining three terms constitute the selection 
free energy: respectively, thermal selection, dilution selection 
term, and finally Fref linear in the effective fields caused by 
dilution. The disorder average of Fref is clearly zero, but 
. such terms give "random exchange fields", so called as they 
(may) act like a random field on the discrete order parameter 
remaining after selection^' 

Next we average Fsc\ over disorder. Recall from ( 12.51 ) that 
h = JP-^e and (eiEj) — Sij{e^) allows us to average AF = 
Fdii [given by dinii]: 
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(Fdn) = --(e^)Tr 



PJfA^ A^P^B^PA MJP^ 



(2.19) 

Substituting J hoi - A from ( 12.101 ). Eq. ( 12.191 ) expands 
into zero, first, and second order terms in Hq. By repeated 
use of ( 12.91 ) and ( 12.14b . the contributions of the two terms in 
parentheses inside ( |2.19t are seen to cancel in both the 0{ho) 
and 0(/i§) parts of the expansion. The result is 



TrfPAP^ B 



(2.20) 



Now Tr(PAP^) = Tr(AwxAr) = Nho, which is the same 
for all ground states. Gathering terms, the averaged selection 
free energy is 



(i^sci) = ifcsrTr(lnB) 



1 



e')Tv(B-')--(e')Nho, 



2' ' ^ '2 

(2.21) 

where we have used ln(dct B) = Tr(lnB). 

It can be shown that ( |2.20t vanishes for collinear states and 
is negative semidefinite (see AppendixlAb. Therefore dilution 
must favor noncollinear states. 
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B. Effective Hamiltonian for Selection 



where f^J""^ is the 0(((5B)") term. In ( IZ27b , 



In this section we re-express the classical selection (free) 
energy ( 12.211 ) in the form of an effective Hamiltonian. In this 
form, it is easy to join the selection terms with other pertur- 
bations that also tend to select states, such as anisotropics or 
magnetoelastic couplings. The selection term takes the bi- 
quadratic form ( 11.2b . for all exchange-coupled frustrated vec- 
tor magnets. 

In the J1-J2 antiferromagnets on the square lattice [(7r,0) 
order], the fee lattice (type II order) or the bcc lattice (bcc type 
II order), the ground states consist of sublattices which can ro- 
tate independently of each other (provided order-by-disorder 
is neglected!). These sublattices are connected only by Ji in- 
teractions, so J1/J2 can be used as a small parameter Indeed, 
the lowest-order term in the J1/J2 expansion does have the 
biquadratic form (11.2b for the J1-J2 XY model on the square 
lattice'^* or for quantum fluctuations in the large-5 Heisenberg 
model on the bcc"^. However, for many other cases (e.g. fee 
type I or type III), the degeneracies entail correlated changes 
in different sublattices, so there is no natural small parameter. 

We deploy a different trick, which has no guaranteed small 
parameter. It is based on the central role played by the matrix 
B in the previous subsection. Recall (12. 21b gave the selection 
free energy in terms of Indct B and TrB^^. The matrix B 
was defined in real space by (12. 14b as 



B^j — Si • SjGi 



where 



(2.22) 



(2.23) 



is a Green's function. (For this section, we drop the super- 
script in S^^-*; it is understood that is a mean direction.) 
Although Gij has the translational invariance of the lattice, it 
is not simply a function of — rj if the sites do not form 
a non-Bravais lattice (e.g. honeycomb, garnet, or pyrochlore 
lattices). Since the A matrix is a discrete approximation of the 
Laplacian, Gij is likely to behave at long distances as 1/ r^j 
(in d = 3) or as In r^j (in d = 2). 

Now for the trick: let B be the average of B over the whole 
ground state ensemble. 



B, 



(S.i ■ S,)cGi 



(2.24) 



where the angle brackets refer to an unweighted average over 
all ground states.^'' (Notice B has the full symmetry of the 
lattice, higher symmetry than B has.) We shall formally con- 
sider 



(5B = B - B 
to be our small parameter Thus 



<5B,; 



Gij [Si ■ Sj 



(S, • 



(2.25) 



(2.26) 



Then from ( 12. 21b by expanding the logarithm of B (for ther- 
mal selection) or its inverse (for dilution selection) in powers 
of (5B, we arrive at 



(2) 



(2.27) 



/sol'"' - ifcBTln(detB) + i(e')Tr(B-i), (2.28) 



f ,(1) - i V r^s- s 

/sol — 2 '^ij ''^ "^3 ' 

i,3 



(2.29) 



where 



(2.30) 



and 



with 



/''^ = -7 E ^y-.«(S, • Sk){Si ■ so, (2.31) 



K 



jk,li 



GjkGii 



kBT{B-%{B-') 



kl 



+ (62)(B-2),,(B-i) 



kl 



(2.32) 



Here "symm" means the expressions are to be symmetrized 
over all distinct ways of pairing indices, such that each pair 
has one from jk and one from li (two or four ways, respec- 
tively, in the two terms). Each term includes contributions 
from both Ft and Fdn. 

Clearly /^"^ does not cause selection, since it produces only 
constants dependent on T and (e^). As for the /(^' term of 
( 12.32b . it takes the form of a renormalization of the ordinary 
exchange couplings Jy — > + SJij, where 



5Jij = 2kBT{B-^] 



ij Gji 



(B 



%G,,. (2.33) 



2J 



(Here SJij has the full lattice symmetry, since B and Gij do.) 
Within the ground state manifold, /^^^ is normally constant 
(assuming the degeneracy is generic by symmetry, and not due 
to a particular ratio of the couplings). Nevertheless, the renor- 
malization will shift the phase boundary between the kind of 
order in question and a competing spin order. 

When can we simplify Z'-^-'? Commonly, and certainly 
when the sites form a Bravais lattice^"'^^ one can write any 
classical ground state as a linear combination of the (degen- 
erate) optimal eigenvectors of the matrix, the so-called 
Luttinger-Tisza construction^^'. The (vector valued) coeffi- 
cients play the role of order parameter components {N^}. 
Obviously, /(■^^ will be quartic in {Nq}. If group theory tells 
us there are m independent quartic combinations - m = 2 for 
the type III fee antiferromagnet'^ - then the parameter space 
of possible f^^^^ functions is only m dimensional. We could 
exactly reproduce these functions using ( |1.2b out to the m-th 
nearest neighbor. 



C. Real-space perturbation approach 

There is an alternative, cruder path to the biquadratic form, 
in the case of classical thermal fluctuations^^. It is included 
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here because it is fully analogous to the approach we use for 
quantum selection in Sec. Hill 

The idea in this approach is to treat the site-diagonal part 
of ( 12.4b as the zero-order part, and to pretend the inter-site 
terms are a small perturbation Ti' . (In fact the site-diagonal 
and inter-site coefficients are manifestly the same order in the 
Ji j's. However, as elaborated in Sec. lIIIBl we may consider 
l/z to be a sort of small parameter, where z is the coordination 
number) 

We continue to assume that the local field /iq is the same 
on every site. Thus our zero-order Hamiltonian is quite trivial 
and tractable, and the perturbation free energy is derived from 
the cumulant expansion to second order: 



(2.34) 



where (...)o means the expectation in the ensemble with just 
hf). (I omit the terms in {H')o, since {H')o — 0.) 

By taking expectations in this trivial ensemble, one 
sidesteps most technical annoyances due to the different ori- 
entations of spins, which necessitated the auxiliary integration 
in ( 12.121) . We simply have 



ha 



(2.35) 



where fij = S^. (This notation for the unit spins is consistent 
with SecHn]) Then 



Terms mixing two different pairs vanish since {(Ti)o 
(|236] |. 



\<^ili<^jfj.O'iiyO'ju)0 



l + (n, -fi.f 



\ ho J 

Putting it all together, we get the form (ll.2l l with 

J? 



J2I 



(2.36) 
= 0. In 

(2.37) 
(2.38) 

(2.39) 
(2.40) 



That is exactly what we would get from ( 12.32b if we took 
Aij = hgSij in place of ( 12.101 ), which would be consistent 
with our pretence that J is small compared to Hq. 



III. HARMONIC QUANTUM FLUCTUATIONS 

A similar effective Hamiltonian emerges from an approxi- 
mate treatment of the selection due to quantum spin fluctua- 
tions within the spin-wave approximation. The basic idea of 



the method was used by Long^ for calculating quantum se- 
lection in type I FCC antiferromagnets. We reformulate this 
approach to be explicitly independent of the "gauge" arbitrari- 
ness in choosing local frames (Sec. Ill B), and show that quan- 
tum fluctuations always favor collinear states. 



A. Set-up for harmonic Hamiltonian 

We start from a classical ground state spin configuration, 
defined by a set of classical directions {ui}, and set up the 
Holstein-Primakoff representation in the standard fashion (for 
noncoUinear spins^*"). In contrast with Sec. [Ill from here on S 
denotes the spin operator and is not normalized to unit length. 
A local orthonormal triad (xj, y^^ z^) such that Zj = fij, and 
we quantize along the Zi axis; we let a bar distinguish spins 
written in this local basis. Holstein-Primakoff bosons are in- 
troduced as usual: 



Siz — S a - ctj, 



2S 



(3.1a) 



(3.1b) 



where Si+ = Sj_ . (We restored the explicit powers of S that 
were omitted in the classical discussion of Sec. |ll]) Of the 
Taylor series for the square roots, giving the semiclassical ex- 
pansion in l/S, we just need the lowest order approximation, 
Si- « \/25aJ. That means 



(e+al 



(3.2) 



Here e± = (x ± iy) / \/2. Next (13.2b must be substituted into 
the spin Hamiltonian dl.lb . The zeroeth order term (in {a|} 
and {fli}) is the classical ground state energy, while the first 
order term vanishes since the ground state is classically stable 
(each hi is lined up with the local field h^). So we only treat 
the second order, i. e. harmonic terms. The resulting quadratic 
Hamiltonian is a quantum analog of (12.41 ): 



harm 



= const — hi(S — itti) 

i 

+ i ^ Jij (e+a] + e-ai^ + e_aj^3.3) 



The usual course is to (Bogoliubov) diagonalize ( 13.31 ), giving 
the harmonic spin wave modes and frequencies {ujm}. Then 
the harmonic zero-point energy is 



£^haim({fii}) = E ^hUn 



(3.4) 



Since the spin-wave spectrum is usually different for 
symmetry- unrelated states, ( 13.4b typically breaks the degen- 
eracy with the energy scale 0{JS) (down by l/S* compared 
to the total classical energy.) As in the previous section, we 
desire simple explicit expressions for the dependence on fi^ 
that is implicit in (13.4b . i.e., an effective Hamiltonian. 

To implement this for a frustrated or highly frustrated anti- 
ferromagnet requires 
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(i) The assumption - not always valid - of a periodic 
state; in any case, the magnetic unit cell is typically 
large leading to large matrices; 

(ii) The choice of an expUcit frame for each site, normal 
to its spin, as the basis for possible spin deviations, with 
its gauge arbitrariness^^; 

(iii) the calculation should be repeated for every one of 
its (continuum of) inequivalent ground states. 

As in the classical selection calculation, all this makes the cal- 
culation cumbersome and sometimes impossible. 



B. Real-space perturbation approach 

To get an approximate answer of simple form, we shall 
decline to diagonalize the Hamiltonian, and instead treat the 
inter-spin terms in (13.31 1 as if they were small perturbations: 
that is, Tio is defined to be the first term, and H' is the second 
term. (This was pointed out by Long**.) 

The ground state |0) of Ho is unique and obviously Siz — S 
for all sites; the first nontrivial term must come from second- 
order perturbation theory. The excited-states in 7i'|0) have 
two spin flips on coupled sites with an excitation cost 

S{hi + hj), as created by the term with two boson creations. 
(In most frustrated systems, hi = ho the same on all sites, 
allowing for simplifications.) Thus, in second order perturba- 
tion theory, we get the correction^ 



6E- 



Jf. 



2 ^ h. 



hi 



-i+ ■ I 



(3.5) 



Now we take advantage of the fact that {ici, y^, z^} are rows 
of a (proper) orthogonal matrix: 

Xi\x^, + yi\yif_, + ZixZii_, = (5a^; (3.6a) 
Xixyt^i, - yixxtfj, = e\fj.uZi,y (3.6b) 



of order S and of course hi = S\ J^j Jij^jl is manifestly of 
order Jy . The implicit parameter of this second approxima- 
tion is evidently "J/ /iq" where J is an appropriately weighted 
average of the couplings, and ho is the averaged local field on 
(if that is not already uniform). Then ho ~ Jz, where z is an 
appropriately weighted count of the nearest neighbors. Thus, 
we are implicitly making large coordination number approx- 
imation. It should work best on the fee lattice (z = 12), or 
in lattices with J2 comparable to Ji. One expects it to be 
poorer quantitatively in (say) the nearest-neighbor Kagome 
lattice {z = 4). However, a biquadratic form with a fitted 
coefficient might still be a reasonable approximation of the 
actual functional form in that case. 

Expanding the square in (13.8b gives a constant term, a renor- 
malization of the coupling due to fluctuation so (in terms of 
the unit spins {fii}) 

^ (3.9) 



Tunit J q2 



Aho 



Jf. 



and most importantly a term of the same biquadratic form as 
introduced for thermal and dilution selection, Eq. (I1.2l i with 



S J ^, 



(3.10) 



Because of this form quantum fluctuations must favor 
collinear states. This conclusion extends to all degenerate vec- 
tor magnets with site-independent ho- 



C. Comparison to alternate approaches 

As mentioned above, in the case of independent sublat- 
tices, one obtains an effective Hamiltonian by perturbation in 
J1/J2, where Ji is the inter-sublattice coupling. In the case 
of the Ji-J — 2 square lattice antiferromagnet, the spin- wave 
dispersion is easy enough to get analytically, namely 



huj{q) = S<[C{q)-2Ji[cosq.^+cosqy]\(3.ll) 



Hence, assuming the summation convention. 



ix\ - iyx){x'^ - iy'x){x^, - iy^,){x'^ - iy'^) 
(1-z-zf (3.7) 



Substituting into (13.5b . we get the main result 



o2 72 



2 hi + hi 

1,3 



(1 - fi, • fij 



(3.8) 



Is there a small parameter in this expansion? The small pa- 
rameter which justifies keeping just the harmonic-order term 
is 1/ S, as always in a Holstein-Primakoff expansion, and so 
the selection term is of order 1/5 relative to the classical en- 
ergy, expansion at harmonic order But we made a second, 
more brutal approximation of the harmonic-order in treating 
the second term in ( 13.3b as a perturbation: both terms there are 



1/2 



x(C(q) — 2J I COB ^[cos qx — cos qy] 



(3.12) 



where C(q) = 4J2(1 ^ cos g^: cos ) . Integrating the zero- 
point energy over the zone and expanding to second order, we 
get the spin wave energy (per unit cell) 

t2 

E.^ = E^l^ - 77^ G-i cos2 $ (3.13) 



8J2 



where 

G-i 



cPq (cos q^ — cos g^)^ 
(27r)2 1 — cos (73; cos (7j, 



0.727. 



(3.14) 



If we match this with (11. 2b (with two bonds per unit cell), 
noting cos^ $ ~ (fi j • rij) for nearest neighbors, and remem- 
bering the local field is /iq/S* = 4J2, we get 



4/10/5 



G_i. 



(3.15) 
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Comparing to ( 13.1 Oi l, we see the real-space hopping approach 
overestimated the biquadratic selection energy by ^ 33%, 
which is not bad in view of the crudeness of the approxima- 
tion (and considering that the coordination number is not so 
large for this lattice.) 

Two other tricks have been used, instead of the real-space 
expansion trick, to arrive analytically at an simple effective 
Hamiltonian of biquadratic form The first trick is expanding 
around an infinite-range modeP''. Comparing to that helps il- 
luminate the derivation in this paper (Sec. 1111 B I above), since 
either approach is a way of setting up mean-field approxima- 
tion. That means breaking up the Hamiltonian as Hq + H' 
such that Tig has a trivial (and tractable) state that approxi- 
mates the ordered state of interest, and expanding in H'. 

Just as in statistical mechanics, there are two standard ways 
to set up the mean-field theory in real space. In the approach 
of Sec. 1111 B1 Ha consisted of local fields that fix each spin in 
its classical direction, and the wavefunction is a direct prod- 
uct of spin coherent states. On the other hand, in Ref. [3^ 
Ho consists of infinite-range spin-spin couplings Jij ; they de- 
pend on the sublattice of i and j, but are unchanged if ei- 
ther site is translated by any number of (magnetic cell) lattice 
vectors. The latter approach, then, is limited to states with a 
known kind of long-range order, such as face-centered cubic 
antiferromagnets, and cannot be used for highly frustrated an- 
tiferromagnets, which have a huge ensemble of non-periodic 
classical ground states. (Strictly speaking, the infinite-range 
approach is not even well posed if this ordered state gets mod- 
ified into a texture with slowly varying directions, or if it has 
thermal spinwave fluctuations. Thus, the advantage of the 
local-fields formulation is that it handles arbitrary configura- 
tions of classical directions {n;}. On the other hand, an ad- 
vantage of the infinite-range formulation^^ is that it explicitly 
obtains finite-size corrections in the spin-wave selection en- 
ergies. (The corrections are big at system sizes tractable by 
exact diagonalization.) In the systems where both approaches 
apply, they give the same answer^^. 

A second trick works in the case of a pyrochlore lattice (or 
any other lattice made of corner sharing triangles and tetrahe- 
dra). The sum in (13.4b can be represented as the trace of the 
square root of the dynamical matrix (basically H'). A Tay- 
lor expansion of that square root gave a series of terms which 
depend on the classical directions of each spini^^i^ 

That was only applied to collinear ground states, but (see 
Sec. V of Ref. 34) it works with small modifications for any 
ground state configuration. It seems that the second order term 
of that expansion is identical to the result we found here. (In 
the collinear case, this term was trivial, and nontrivial effects 
were found only in higher-order terms representing loops in 
the lattice.) The trace viewpoint, however, is more powerful 
in that it allows computation of higher terms; also, in some 
cases partial resummations of higher terms might give better 
numerical approximations for the coefficient in the effective 
Hamiltonian, analogous to the resummations in Ref.[34l 



D. Kagome and similar antiferromagnets 

For some frustrated magnets, the effective Hamiltonian 
(13.8b takes the same value in every classical ground state and 
hence does not break the degeneracy. In particular, this hap- 
pens for lattices (such as the kagome) built from triangles, 
each of which has a 120° spin arrangement in a ground state. 
Since hi hj = — 1/2 for all nearest neighbor pairs, Eq. ( 13.8b 
reduces to a constant. 

In such systems, coplanar states are the closest thing to 
collinear states within the ground state manifold, and these 
indeed have the lowest harmonic zero-point energy. 

What would be the natural terms in an effective Hamilto- 
nian, analogous to ( 13. 8b and ( 11. 2b ? Since three spins are re- 
quired to distinguish a coplanar from a noncoplanar state, one 
might think a term is needed of form 

|n, • hj X fifep = 1 - (n, • hjf - (fij • fife)^ - (fife • hif 
+ 2{h,-hj){hj-hk){hk-hi). (3.16) 

The simplest spin term which selects for the coplanar states 
is just a second-neighbor biquadratic term of form Eq. (11.21 ). 
(In a nearest-neighbor Heisenberg model, that would not get 
generated by ( 13.8b .) On the kagome lattice (and other lattices 
of corner-sharing triangles), the third-neighbor coupling 
(for spins in the same row, not across a hexagon) must be equal 
to the second-neighbor term K2, for two reasons: 

(i) if it were not, this would split the degeneracy be- 
tween different coplanar states, whereas we know they 
are exactly degenerate according to the harmonic-order 
selection energy 

(ii) at the perturbation order that this term mostly comes 
from, there is no distinction between paths to these two 
neighbor-of-neighbor sites. 

An extension of the Sec. 1111 B I calculation to fourth order 
perturbation theory does generate the desired neighbor-of- 
neighbor biquadratic terms on the kagome lattice. (The trace 
form from Sec. V of Ref. [13 is probably an easier way to set 
this up.) Specifically, K2 and come from processes in- 
volving a string of sites i, j, and k (not on the same triangle). 
There is a spin exchange involving the pair (ij), another in- 
volving (jfc), and then the same two pairs (in either order) so 
that we return to the ground state of the trivial Hamiltonian 
H^. 

The phenomenological term adopted in Ref. 38 is the same 
term that would follow from the second-neighbor biquadratic 
interactions. (This was assumed as the barrier potential to esti- 
mate a spin-tunneling amplitude in the kagome antiferromag- 
neta) 

Unfortunately, this functional form is wrong for the 
Kagome case^'', and for coplanar states in general. Instead, 
the effective Hamiltonian scales non-analytically as |sin0|, 
where 9 is the angle between the spin planes in adjoining 
parts of the lattice^'*. Sec. V of Ref. |34| has clarified why 
the spin-wave energy cost, when spins rotate out of one of 
the discretely selected ground states, scales linear in angle 
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deviations for the (coplanar) Kagome case but quadratic for 
the (colHnear) pyrochlore case. This cusp should get some- 
what rounded once anharmonic fluctuations are taken into ac- 
count^**. It is an interesting, unsolved challenge how to for- 
mulate the calculation of a local effective Hamiltonian, in the 
spirit of (13.81 1. that captures the cusp behavior 

Nevertheless, a phenomenological form favoring copla- 
narity, such as a neighbors-of-neighbors biquadratic term, 
would be an improvement over classical simulations that do 
not account at all for quantum effects. In static correlations, as 
can be accessed by Monte Carle"**', one expected effect would 
be a much more robust stabilization of coplanar states (at 
much higher temperatures). In the dynamics, as can be simu- 
lated by molecular dynamics'", the coplanarity term gives a fi- 
nite frequency to the band of zero-frequency spin wave modes 
(except the Goldstone mode). 



IV. DISCUSSION 

To summarize, we have derived effective Hamiltonian 
terms for selection free energies from spin wave fluctuations 
due to thermal excitation, quantum zero-point motion, or to 
dilution disorder; the result has the form of an effective bi- 
quadratic Hamiltonian ( ll.2l l in every case, with the coefficient 
proportional to temperature, 1/5, or dilution concentration, 
respectively. In each of the three cases, our key trick was to 
rework the problem until the spin configuration became the 
entries in a matrix, and then expanding in them. The same 
idea has, in one guise or another, been successful for several 
independent examples of getting an effective Hamiltonian by 
integrating out fluctuations. Besides Ref. 2, which is a precur- 
sor (for the classical undiluted fee lattice) of our calculation 
in Secini this was done for quantum'*^ or classical'*^ kagome 
systems, as well as pyrochlore antiferromagnets^'*'*'*. 

Sometimes, it is possible to compute numerical energies for 
a large database of configurations, and then fit an effective 
Hamiltonian empirically'*"*. Even that brute-force approach 
depends critically on analytics which suggest the proper func- 
tional form to be used for fitting. 

Now we turn to the applications of effective Hamiltonians. 
They are convenient for predicting the ground state spin pat- 
tern, when the Hamiltonian is complicated by too many kinds 
of interactions, e.g. spin anisotropics of all sorts, dipolar in- 
teractions, distant-neighbor couplings, external fields, or mag- 
netoelastic coupUngs; also, defects, inhomogeneities, bound- 
aries and domain wafls. Thus, they allow quick understanding 
of phase diagrams which may have many parameters. For ex- 
ample, phenomenological biquadratic term representing the 
quantum fluctuations was used^ to explain a plateau in the 
dependence on external field of the incommensurate order- 
ing wavevector of CsCuCls. Another one-dimensional model 
with a biquadratic term was studied in Ref*^. 

Effective Hamiltonians have also been combined with 
semiclassical tunneling theory to understand tunneling pro- 
cesses, when degenerate- ^ ^•'*^'"*^ or highly degenerate^ ground 
states are separated by barriers. 



A. Pitfalls of classical simulations 



As one of us noted in Ref. [Tsl large S justifies visualiz- 
ing each spin as a fixed-length vector, but it does nof justify 
a purely classical simulation of the system, as is commonly 
don o^'''*°i'**'^'''^'i^^'^^ . Notwithstanding the well-defined spin 
directions, the unmodified classical Heisenberg model gives a 
qualitatively wrong picture of the behavior. 

The reason can be expressed in terms of different energy 
scales, as outlined in Ref. (49t The highest scale is Ej ^ 
JzS^, the scale of the mean-field ordering temperature and 
Curie-Weiss constant; when T ^ Ej it is already a good 
approximation that the system is in a classical ground state 
with small fluctuations. Another scale is the temperature T* 
where ordering, freezing, or other phase transitions take place: 
by definition, a "highly frustrated" magnet' is one in which 
r*/TMF<i; the ratio can approach 10~^. We are usually in- 
terested in lower temperatures T ^ T*, where the interesting 
changes occur 

The scale which has been ignored is that of harmonic spin- 
wave fluctuations, Ecow ^ JS < Ej. As long as Ecou < 
T < Ej, the system is certainly classical. But we're most 
interested in T T* ^ Ecou, where a classical descrip- 
tion is not valid. Consider for example'^ the energy barriers 
against flipping from one collinear (or coplanar) state to an- 
other: those due to spin-wave zero-point energy are greater 
(by the ratio Ecow/T) than the classical barriers. That will 
drastically affect any thermally activated dynamics, and will 
enhance the ordering tendency (e.g. a large S pyrochlore anti- 
ferromagnet will develop a long range order at some T > 0.). 
Indeed, in the T < Ecow regime where the classical fluctua- 
tions are just a small correction to the quantum fluctuations, 
the T-dependent selection free energy does not have the form 
it does in a purely classical modeP'*. 

I believe there is an easy fix: if the phenonmenological term 
jl.H (or whatever is appropriate for the system in question) is 
added to the classical Hamiltonian by hand, the modified clas- 
sical simulation can closely represent the low-temperature be- 
havior Of course, such a simulation are much tractable (in 
terms of temperatures, system sizes, and observables) than 
quantum Monte Carlo This could be incorporated both into 
Metropolis Monte Carlo (MC), for modeling the thermody- 
namics and phase transitions, as well as molecular dynamics, 
for modeling spin-wave excitations or magnetic relaxation be- 
haviors. (If quantitative results are needed, it would be best to 
empirically fit the coefficients in ( 11.21 ) first.) Occasionally a 
more exact selection free energy has been used in MC simu- 
lations^^, but that is enormously cumbersome (requires diag- 
onalizing an 0{N) dimensional dynamical matrix every MC 
step.) Having the form of a local effective Hamiltonian makes 
it very easy to allow for the quantum effects, which is neces- 
sary if the simulation means to represent a Heisenberg magnet 
realistically. 



9 



Acknowledgments 

This work was supported by NSF grant DMR-DMR- 
0552461 and its predecessors. We thank Assa Auerbach, Ev- 
genii Shender, Uzi Hizi, Roderich Meissner, and Tom Kaplan 
for helpful conversations. The original work was completed 
at Boston University. 



APPENDIX A: ALTERNATIVE METHOD OF DILUTION 
SELECTION 

We explore an alternative method of calculating the dilution 
selection energy E^-a, using Lagrange multipliers to enforce 
the unit length constraints. Rather than the model dilution 
employed in Sec|II] we consider the more realistic case 



(Al) 



where rjij = e^ej — for [i ^ j). Here = 1 with probabil- 
ity p and is otherwise zero. Thus {rjij) — 0. Then, 



where 5H = \ Y^i.-, ^v^^i ■ (^j^ and 



(A2) 



SH' = E„- -Nhop{l -p) + SHnEF + SHi. (A3) 



The random exchange field contribution is 



1 



2^- 



(0) sj(0) 

3 ' 



(A4) 



which averages to zero and is ignored as before for computing 
selection, and 



Here 



SHi = ^cTi • hi. 



hi — ^ ^ rjij S^j ^ , 



for realistic dilution, or 



(A5) 



(A6) 



(A7) 



for the model dilution of i2.2\ . To minimize the total zero 
temperature energy with respect to the spin deviations <Ti for 
a fixed configuration of the quenched disorder variables rjij, 
we need to include a Lagrange multiplier Ai to ensure that the 



unit length constraints cri ■ s|"^ = are obeyed, yielding 



(0)n 



(A8) 



Introducing a matrix notation as in Sec|ll] the Ai are deter- 
mined by the constraint equation to be 

A = -B^PA^h. 



When (|A9l l is substituted into (IA8l l. and the result substituted 
into (IA2b . we obtain an expression for the total energy in the 
presence of the specific configuration of dilution variables: 



1 



p--^(h^ A-ip^B-iPA-ih - h^A-ih)(A10) 



identical to (12.171 1. 

Both terms take the form h^Mh, which must now be aver- 
aged over the disorder. In components the average becomes 

J2 E M^,,,uSi^^sj^JAaA,k{mrlJk)■ (All) 

fJ-f ijkl 

The average is 



+ p'^{l ~ pfiStjSik + SikSij), 



(A12) 



if we also replace the A's in ( lAllI ) by J's to ensure that i ^ 
l,j k. The terms proportional to (1 — p) represent effects 
of independent missing spins (defects), and will be seen to 
reproduce the form of i?dii calculated in Sec|II] The (1 — 
pY terms account for pairwise correlations between different 
defects. Thus, 

(h^Mh) = p'''(l-j3)Tr[/ioPJMP^ + /ioPMJP^ + 
PJMJP^ + /igPMP^] +p^{l- pf X 
X ^ J^,J^,s'i} [sf^M,^,,, + (A13) 

It is now straightforward to evaluate i?dii by substituting 
(IA13b in ( lAlOl i. Note that the second term in (IA13l l is smaller 
than the first by order 1/ z where z is the coordination number 
(z = 12 or 18 for face-centered cubic, depending whether 
J2 is included.) After performing simplifications we find the 
dilution energy 

1 



E^,, = -p{1-p)[TtB-' ~Nho] 



+ 1(1 ~P? 



JijJijS, 



(0) 



S. 



(0) 



[A-ip^R-ipA-i -A-^] . . 
c(o) [A-ip^B-ipA-i-A-i]. 



(A14) 



(A15) 



(Some simplifications of the 0(1 —p)^ terms are possible, but 
the result cannot be written in terms of conventional matrix 
operations.) 

Clearly the first term in ( IA15I ) is just the dilution part of 
( I2.2II 1 of SecHIlif {e^) is identified with (after rescal- 

ing Jij Jlj = JijP^). As ai-gued in Sec|ll]i?dii — for any 
collinear state. Here this is clearest from examination of ( IA8I ). 
In a collinear state {T — 0), all the S^"* may be taken along 
(A9) the z axis. Then clearly the only allowed directions for ai 
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fluctuations are in the xy-plane, eliminating the second two 
terms from ( |A8t and forcing cTj = for all i. Clearly then 
i?dii = 0. On the other hand the form of p^5H + SHi shows 
that in a generic noncollinear case, where cTi ^ 0, the dilu- 
tion energy for every realization of the quenched fluctuations 
will be reduced — the system takes advantage of the fluctu- 
ation. Thus for a generic noncollinear state, we must have 
Edii < stricdy. 

Systems with an exceptional degeneracy may have i?dii = 
even for noncollinear states. For example, the 2D XY square 
lattice AF with J2 < and Ji = 0. Then the system forms 
two decoupled simple antiferromagnetic sublattices. The cou- 
pling from the dilution then comes only from spins on the 
same sublattice, which are collinear Thus i?dii = even 
though the total spin configuration is noncollinear 

It should be emphasized that the terms above are not the 
only contributions of 0(1 — p)^, since neglected terms of or- 
der cr^ also contribute. Only in a case where some small pa- 
rameter controls the size of the cTi will the exhibited terms 
give the dominant 0(1 — p)^ contribution. 

APPENDIX B: EFFECTS OF MAGNETOELASTIC 
COUPLING 

It is well-known that exchange striction, an isotropic effect 
of nonzero magneto-elastic coupling, also leads to an effective 
biquadratic exchange. It is long known these favor collinear 
states, e.g. in MnO, and they play important roles in highly 
frustated spinels^*. To see how the biquadratic form arises in 
this case, we add to the Hamiltonian the elastic energy H^i 
and magneto-elastic coupling i/mo- The elastic energy is 

-ffoi = ^ f (fixe^xEf,,, . (Bl) 



Here e^A is the strain tensor, and {cka^j,} are the elastic con- 
stants, separately symmetric in (kA) and (/ii^), and under their 
interchange. (We use the summation convention for compo- 
nent indices.) The magneto-elastic coupling is given by 

H^e^lY. eKA«:Aa,^^S, • S, , (B2) 

obtained simply by expanding the exchange constants in lat- 
tice displacements. Thus the magneto-elastic constants are 
defined by 

< = -i^^^).7^ (B3) 

Next the energy Hc\ + i/mc is minimized with respect to the 
strains Ski and the result substituted back in this expression. 
The result is 

Hgj = EE[^"']«v^"'^<r(s. • s,)(Sfc • si) . 

ij kl 

(B4) 

This has the biquadratic form. Within the ground state man- 
ifold (IB4l i has an expression in terms of the same quartic in- 
variants mentioned at the end of Sec. IllBI 

To estimate the size of this effect, we assume the exchange 
constants decay rapidly with distance, according to J(r) = 
Joir/royf,^'^ yielding 

«-(^J(r) . (B5) 

There is some evidence that for superexchange, / is between 
10^ and 14,2 
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